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Abstract 

. General ideas of gauge/gravity duality allow for the possibility of time dependent 

solutions that interpolate between a perturbative gauge theory phase and a weakly 
curved string/gravity phase. Such a scenario applied to cosmology would exhibit a non- 
geometric phase before the big bang. We investigate a toy model for such a cosmology, 
whose endpoint is the classical limit of the two-dimensional non-critical string. We 
discuss the basic dynamics of this model, in particular how it evolves toward the double 
scaling limit required for stringy dynamics. We further comment on the physics that 
D , will determine the fluctuation spectrum of the scalar tachyon. Finally, we discuss 
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various features of this model, and what relevance they might have for a more realistic, 
higher dimensional scenario. 

1 Introduction 

The principle of holography P] 121 IH] seems to imply that space-time geometry emerges 
dynamically from more fundamental degrees of freedom. In the AdS / CFT correspondence jU 
E], for example, it is the CFT that enjoys a more fundamental, non-perturbative description, 
while the stringy, geometric description is applicable only in the regime of large 't Hooft 
coupling. For time dependent solutions, it is possible that the geometric description is only 
applicable for part of the time. In particular, one may consider a cosmology that will be 
qualitatively very different from existing proposals, where the universe at early times does 
not enjoy a geometrical description. 
*novak@ihes.fr 
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Using the qualitative features of the AdS/CFT correspondence as a guide, we assume 
that the correct description of quantum gravity has a holographic description given by the 
dynamics of some theory of N x N matrices. This theory may be either a regular large 
iV gauge theory or something more general, such as a suitable open string field theory 
with iV-dimensional Chan-Paton factors. We further assume that in the large iV limit this 
theory has a (closed) string expansion with an effective 't Hooft coupling A e ff oc N. 
When the 't Hooft coupling is small the theory enjoys a (suitably generalized) perturbative 
gauge description, while large 't Hooft coupling corresponds to a weakly curved gravitational 
description, similarly to the case in AdS/CFT where the AdS radius R is proportional to A^ 4 . 
Now if we consider an expanding universe, suitably chosen holographic screens should also 
be expanding. We take a simple interpretation of holography, wherein if a comoving region 
of spacetime expands, the growing entropy limit implies a growing number of fundamental 
degrees of freedom describing the region. Thus we are imposing that locality continues to hold 
in some generalized sense. As the universe cools and expands, an increase in the number 
of degrees of freedom implies a growth in the number N characterizing the fundamental 
"gauge" description. Thus we assume that out of some N tot x N tot total degrees of freedom 
(with N tot perhaps infinite), and at a particular epoch t, all but a certain N e g(t) x N e g(t) 
subset of the degrees of freedom are frozen out from the dynamics. The cooling of the 
universe then involves the thawing out of these other degrees of freedom, which then allow 
for the holographic description of the growing universe. Since we are now, in effect, allowing 
for a time dependent N, we expect that the dynamics of the large N expansion is controlled 
by an effective 't Hooft coupling X e s(t) oc N e g(t). As the universe cools, from a small subset 
of highly excited degrees of freedom to a much larger set in the future, the 't Hooft coupling 
grows, and we travel from a weak gauge coupling regime to a strong coupling regime. It is 
the strong 't Hooft coupling regime which should have a stringy, geometric description. 

As a potential mechanism for cosmology, we merely point out some basic ideas, leaving 
a more detailed investigation to future work [7j . 

• We begin with a very high temperature, weakly coupled gauge theory. This system 
will have some characteristic spectrum of fluctuations. Further, depending on the order 
of the transition between the gauge and gravity phases, there will be some character- 
istic spectrum. For example, in condensed matter physics, second order transitions 
produce scale-invariant fluctuations in effective field theories. These fluctuations will 
then naturally translate into some set of fluctuations in the gravity description, and 
in particular may seed the formation of large scale structure. It is interesting in this 
respect to note that the cosmic microwave background fluctuations are indeed very 
nearly scale invariant, a feature conventionally associated with slow-roll inflation. 
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• If the radius of curvature is indeed proportional to some positive power of the time 
dependent, effective 't Hooft coupling, as in AdS/CFT, then a mechanism like ours, 
driving the system towards large N, will move towards a flat universe. In addition, 
since the evolution of the universe will naturally extend past the big bang into the 
gauge phase, the horizon problem can be naturally resolved. The scenario may thus 
provide an alternative to inflation. 

• Even in the case that there is still an inflationary epoch, our scenario may be useful for 
understanding the initial conditions. That is to say that initial conditions that may 
be very natural for a highly excited, perturbative gauge theory may look rather exotic 
by the time that they are propagated (through the dynamics of the phase transition) 
into initial conditions for the geometric phase. 

• String cosmology is typically studied using the low energy effective field theory, with 
perhaps some exotic non-perturbative features, like D-branes, added to the mix. In 
this context, it has become important to understand the huge space of (apparently) 
allowed string vacua, which has been termed the string "Landscape" [HI E] • in order 
to make solid arguments regarding the statistics of the landscape (as in the work of 
Douglas PUj), it seems necessary to also understand the behavior of strongly coupled 
models (related to the gauge phase), and the motion between strongly and weakly 
coupled regions. 

Most modestly, then, this scenario is important for understanding the pre-inflationary 
phase in Landscape models which might differ, for example, from those that exhibit eternal 
inflation. Perhaps a bit more ambitiously, it may also serve as an alternative to inflation 
and as an explanation for the big bang and large scale structure of the universe. We leave 
this possibility to future work. 

In this paper we make a first step towards exploring this idea within the context of the 
simplest gauge/gravity duality, that of the c = 1 matrix model jTTj. Our aim in this work is 
to see how we would implement our scenario in this simplified set-up, and to see what sort of 
surprises may be in store for us. In section two we outline the basic dynamics of our model. 
In particular, we find that the ordinary c = 1 matrix model does not allow for a suitable 
freezing of the degrees of freedom, and we are led to introduce a mechanism where quantum 
effects (and thus eigenvalue repulsion) are turned on as the system cools. We further find 
that the instability of the bosonic string is crucial in allowing the system to dynamically 
approach the double scaling limit. It is not clear at this point whether these modifications 
have importance for higher dimensional scenarios, or whether they are necessitated by the 
special kinematics of 2D. In section three we study the nature of the gauge/gravity phase 
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transition in our model, and briefly discuss how to extract the fluctuation spectrum for the 
tachyon. We also discuss the most basic features of the two dimensional cosmology. We 
conclude in section four by comparing the two dimensional model and it's cosmology to 
what we may expect to hold in higher dimensions, pointing the way to questions that we 
shall address in future work. 

2 The Basic Scenario 

We begin this section with a brief review of the key features of the c = 1 matrix model 
[TT] . emphasizing it's role as a simple example of gauge/gravity duality, as clarified in recent 
work spearheaded by McGreevy and Verlinde [E] . Since there have been very many papers 
on this subject recently, our treatment will be brief. The model is described by the U(N) 
quantum mechanics of a single N x N matrix M, moving in a U(N) invariant potential 
V(M). In this paper we will be concerned with the bosonic string, which is related to the 
matrix model with a cubic potential. Following a standard recipe, one can integrate out the 
"off-diagonal" components of M, as well as the gauge field, finally reducing the dynamics 
to those of the matrix eigenvalues, which behave as free fermions A» moving independently 
in the potential V(Ai). The ground state of this theory is obtained by simply filling up the 
fermi sea to some fermi energy ep. 

The string theory phase of this system is approached by taking the double scaling limit. 
One fills the fermi sea up to a fermi energy that approaches a critical value e c , the local 
maximum of the cubic potential. Introducing the parameter /i = N(e c — ep), the double 
scaling limit involves taking the rank of the matrix N — > oo and e c — ep — > with fi fixed. In 
this limit the fermion sea near the local maximum describes a two dimensional linear dilaton 
string background. The scalar "tachyon" (which is actually massless in two dimensions) is 
then related to the collective field describing small fluctuations of the Fermi sea. 

Some time dependent solutions of this model, proposed as solutions for two-dimensional 
cosmology, were discussed by Karczmarek and Strominger [X3| - In these models, an incom- 
ing sea of fermions is chosen so that as the fermions move towards the maximum of the 
potential, they form a smooth region of spacetime for some finite time period. This picture 
is qualitatively very different from what we want to produce, as the double scaling limit is 
implicitly taken in the initial conditions. Namely one needs to choose a large enough blob 
of eigenvalues (N large), with momentum tuned precisely so that they get sufficiently close 
to the critical point of the potential. 

What we want to have in our case is a system that dynamically approaches the double 
scaling limit, as the total energy spreads from an initial small subset of highly excited degrees 
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of freedom. If we look at the ground state of the quantum mechanical model with some fixed 
large N, the system already has a particular Fermi energy, and in particular the choice of 
N and €f is such that the system is either close to the double scaling limit or not. Thus we 
need to modify the model so that the quantum mechanical behavior, and thus the eigenvalue 
repulsion, is turned on as the energy spreads among eigenvalues. The dynamics must be such 
that the quantum mechanical behavior is an emergent phenomenon. 

To see how to implement this, we start with the following Lagrangian, where for the 
moment we consider classical degrees of freedom: 



L = Tr -M + -M + -M 6 , (1) 
V 2 2 3 / 

where M is a hermitian N x N matrix. We diagonalize this matrix by writing M = f^Afi, 
where A is the diagonal matrix of eigenvalues Aj and Q is a unitary matrix. The Lagrangian 
now becomes, after introducing the hermitian matrix A = — iClQ^, 



N /1 1 \ 1 N 

^ = E(^ + iA ? + |A3) + I E |4f 

i=l v 7 i<j 



(A, -A,) 2 . (2) 



Finally, it will be more convenient to pass to the Hamiltonian description. Introducing the 
momenta pi conjugate to Aj and conjugate to Aij, we find 



N , s AT , |2 

1=1 x ' i<j v J ' 

Clearly it is the last term that controls the eigenvalue repulsion. Equally clear, however, 
is that with this Hamiltonian the IT-, momenta are constants of the motion. Thus we will 
have to modify the classical dynamics if we wish to have the eigenvalue repulsion turn on 
dynamically. 

More concretely, let us imagine that in our initial configuration, all but one of the eigen- 
values are collected very close together in the metastable minimum at \ = —1/g, and that 
we have one eigenvalue at rest at Ai — N 1 / 3 . The initial eigenvalue displacement is chosen 
so that the initial energy is larger than N/g 2 , which is the order of magnitude of the energy 
of the final, filled Fermi sea. For the bulk of the eigenvalues to be able to stay near the 
bottom, the momenta IT^ will have to begin very close to zero. 

Note that the force between two eigenvalues falls off very rapidly, with the inverse third 
power of their separation. To understand the eigenvalue dynamics with fixed ILj more 
clearly, let's look for a moment at two eigenvalues in the absence of an external potential. 



5 



If we choose the initial conditions (at large negative t) so that the left eigenvalue is moving 
to the right with initial velocity v, lim t ^_ 00 A/ = vt, while the right eigenvalue is initially at 
rest at zero, lim^-oo A r = 0, the equations of motion are easily solved to give 




where II is the "off-diagonal" momentum controlling the repulsion, as above. Thus we see 
that as t becomes large and positive, the left eigenvalue is now at rest at A^ = 0, while the 
right eigenvalue is moving to the right with velocity v. Up to exchange of the eigenvalues, it 
looks simply like the free propagation of two eigenvalues through each other, apart from a 
short period of time when the free eigenvalue separation \vt\ is small compared to U/v. Note 
for future reference that the distance of closest approach for the two eigenvalues is simply 
given by 211/?;. 

With the above Hamiltonian, the initially excited eigenvalue would fall toward the rest 
of the eigenvalues, ultimately transferring the eigenvalue momentum to the rightmost eigen- 
value, which would subsequently escape into the instable region at positive A. All of this 
would happen on a time scale of order 1/yN. We see that in addition to needing the 11^ 
to vary, they need to be able to absorb some of the energy, and thus thermalize the system, 
on a similar time scale. Further, we require that this new thermalization mechanism only 
acts on "distance" scales Aj — Xj < 5, for some fixed S of order 1/Ng. This is required both 
so that we can initially separate one eigenvalue that contains all the energy, and also so 
that we can end up with something like free propagation for eigenvalue perturbations after 
equilibrium is reached and the eigenvalues fill out the fermi sea. To see this, note that once 
the Iljj have grown somewhat, and very high momentum eigenvalues have escaped to the 
region of positive A, that the distances of closest approach, of order Tlij/v will have some 
lower bound, which will ultimately be greater than or equal to S. 

We would like to comment at this point that 't Hooft has argued [14J that reconciling 
holography with ideas of locality suggest a picture similar to what we have here. Namely, 
we find a system with classical, deterministic degrees of freedom, which has a short dis- 
tance effect that both loses information (here through thermalization) and ultimately pro- 
duces quantum mechanical effects (here corresponding to the growth and stabilization of the 
eigenvalue repulsion). We shall discuss this more in the final section. 
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3 From Eigenvalues to Spacetime 



We begin this section by summarizing the picture so far. We begin with all the ITjj very 
close to zero, so that all but one of the eigenvalues can sit at the bottom of the metastable 
minimum, at \ = — 1/g. In particular, we require that the initial spread between eigenvalues 
is much smaller than the special distance 5 introduced above, although the total spread, since 
iV is assumed very large, may be larger than 5. The final eigenvalue begins at rest at an initial 
negative displacement, large on the scale of A 1 / 3 , so that it's energy is large on the scale of 
N/g 2 . The eigenvalue falls down, and as it approaches within the distance S of the leftmost 
eigenvalues, a subset of the eigenvalues starts to thermalize, spreading out the energy by 
increasing some of the Uij momenta. Although the dynamics of our model ensures that 
the eigenvalues remain in the same order throughout, up to eigenvalue permutation we can 
follow the trajectory of our special excited eigenvalue. As it travels through the sea of other 
eigenvalues, it keeps shedding energy through thermalization. Since we imagine that the 
thermalization occurs while eigenvalues are within a distance 5, the process of thermalization 
will have some non-trivial behavior. In particular, initially the excited eigenvalue will have 
a higher velocity, and so will spend a smaller amount of time very close to the leftmost 
eigenvalues, while it will spend more time very close to the rightmost eigenvalues. Therefore, 
in the initial stages of excitation, there will be more energy deposited on the right-hand part 
of the spectrum. Whether this tilt persists all the way into equilibrium will clearly depend 
on how quickly it is reached. 

Although so far our system is still very far from equilibrium, to trace the behavior further, 
it is helpful to look at the equilibrium dynamics of the system, where we now allow the 
momenta to freely vary. The partition function for this system is 



This partition function is straightforward to evaluate. The integral over pi contributes an 
iV dependent constant times (3~ N ^ 2 , while the integral over A {j simply contributes an N 
dependent volume factor. The integral over provides a factor of p- N ( N - 1 )/ 2 ^ anc l a l so has 
the effect of introducing the square of the Vandermonde determinant, A(A) = nf<j(^« ~~ -\?)> 
into the measure for the remaining Aj. Our partition function has thus become 



Z 




i=l i<j \_i=l v 7 i<j 




(7) 
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with Cn an N dependent constant. The integral in the above expression is simply the 
path integral for the zero- dimensional (i.e. no time coordinate) quantum matrix model. 
Now we find ourselves in very familiar territory. Following the classic treatment of Brezin, 
Itzykson, Parisi, and Zuber [T^], we can estimate this integral, in the large N limit, by 
doing a saddle point approximation. The result, after translating to our conventions, gives 
both an expression for the free energy (whose specific form will not concern us here) and an 
equilibrium distribution for the eigenvalues given by 

p(A) = ^( ff -s A )VMM , ( 8 ) 

where the turning points are located at 




A ± = — ±2,l WTWj , (9) 



and a is a solution of the equation 

2g 2 N 



t7(l + tr)(l + 2a) = . (10) 



Note that for our choice of potential, the local maximum is located at A = 0, with the 
metastable minimum at negative A and the instable region at positive A. For our model, 
the eigenvalues will also have a gaussian distribution of momenta superposed over the above 
distribution. If we ignore for the moment the complications of non-equilibrium dynamics, 
we note the following feature. As more energy is deposited into the eigenvalue sea, it's 
temperature rises, and thus j3 will be decreasing. There is a critical point in the behavior at 
(3 = l2^/3g 2 N. At this point A + = a/g and in particular the right edge of the distribution 
will now go as (A + — A) 3 / 2 . We expect this critical point to be related to the double scaling 
limit in our model. In particular, the critical point in the usual matrix model is precisely 
the point at which you expect to get string dynamics, and is also the point when you expect 
eigenvalues to start spilling over towards the instability. 

Returning the discussion to the non-equilibrium dynamics, we will assume that the energy 
transferred to the sea by the excited eigenvalue is greater than that needed to reach the 
critical point. The effect of the dynamics will now be to have a tilt, as discussed above, 
which will affect the momentum distribution of the eigenvalues. Further, we expect that any 
eigenvalues with large enough momenta will escape to positive A, so that as equilibrium is 
approached we have some smaller number of eigenvalues present. 

As long as the friction induced by the short distance thermalization is great enough, the 
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sea will end by filling up the metastable dip, with excess eigenvalues draining off. Here we 
notice another interesting feature of our toy model. In the usual quantum mechanical theory, 
the instability of the bosonic string is a problem, and more precise statements require moving 
to the Type OB string, which has a quartic potential and is thus stable. For our scenario, the 
world-sheet supersymmetric string is of no help, as having a quartic potential would require 
an exquisite fine tuning of the initial conditions. The instability of the bosonic string is the 
precise feature that allows our system to dynamically approach the continuum limit. 

Let us also note the following: the higher the effective friction acting on the excited 
eigenvalue, the more uniformly we might expect the eigenvalue "soup" to heat up. If it 
heats up more uniformly, there should be smaller fluctuations in the momenta and thus less 
"sloshing" of eigenvalues over the edge. This will in turn result in a smaller \x for the final 
state. Noting further that \x is proportional to one over the string coupling (/i oc l/g s ), we 
find that strong friction is correlated with strong string coupling. 

Now we would like to comment on the cosmology seen by a string observer, after the 
system is already sufficiently close to the double scaling limit. In the final stages of ther- 
malization, the sea will still be expanding somewhat. Thus an observer will see some degree 
of spacetime expansion. Since this expansion is fueled by the thermalization, and not by 
normal "quantum" dynamics, an observer may postulate that during this period there is an 
exotic source of stress-energy, which however will be absent at late times. In addition, the 
tachyon, which is given by fluctuations of the fermi surface, will have some characteristic 
spectrum of fluctuations. The details of this spectrum will depend on the initial energy, the 
strength of the friction caused by thermalization, the behavior of the phase transition, and 
the subsequent draining of highly excited tachyons. However, we may note here that the 
phase transition may produce a scale invariant spectrum, while the non-equilibrium dynam- 
ics may superpose a tilt, as sketched previously. We may further expect that very violent 
fluctuations will result in eigenvalues that separate from the collective sea. From the point 
of view of the string theory background, these highly excited eigenvalues should correspond 
to unstable D-branes. In principle, this system is simple enough that this may be studied 
numerically. However, it is not clear at this point if having a more precise answer would help 
our understanding of the scenario as applied in higher dimensions. 

4 Discussion: Lessons for Higher Dimensions 

We will conclude by summarizing what we've learned, and what import it may have for a 
more realistic, higher-dimensional model. Most importantly, we have seen that it is possible 
to implement our scenario in this simplified context; an engine can be set up, operating be- 
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tween a small number (here one) of highly excited degrees of freedom and a large number of 
frozen degrees of freedom, which drives the system to expand into a stringy phase. We have 
also been able to extract some insight into the dynamics that should determine the tachyon 
fluctuation spectrum; this should help in understanding the production of density perturba- 
tions in a realistic scenario. In particular, we see that there are two important features which 
are natural to understand first. The first is what the order of the phase transition is expected 
to be, and in particular whether such a transition produces a scale invariant spectrum. The 
second feature to understand is in which direction the non-equilibrium dynamics will tilt the 
spectra. If these two properties can be robustly understood, then we may be able to falsify 
the picture through comparison with the cosmic microwave background. 

One feature of our toy model that seems rather surprising is that we have had to alter the 
model so that quantum mechanics is an effective phenomenon that emerges at late times. 
It is possible that this modification is made necessary by the special kinematics of 2D. In 
particular, since the extent of the fermi sea is determined simply by fermionic eigenvalue 
repulsion, the only way for the "universe" to expand is to have the eigenvalue repulsion, and 
thus h, growing in time. It is, however, entirely possible that such a modification must hold 
for a higher dimensional scenario, 't Hooft has argued [T3] that reconciling conventional ideas 
of locality with the apparently non-local behavior for quantum gravity implied by holography 
may be resolved in a similar manner. In particular, he imagines a classical theory for gravity 
that is conventionally local. This classical theory, however, has some dissipative dynamics on 
short scales which lead to information loss. Quantum mechanics, and the reduced entropy 
implied by holography, arise through the statistical mechanical treatment of equivalence 
classes of states, where two states are in the same class if they both evolve to a common 
state in the future. Our toy model displays a similar behavior: information is lost through 
the short distance thermalization that allows the off-diagonal momenta to vary. We find 
this prospect very exciting. If the higher dimensional model does turn out to require such 
modification, then it may be fruitful to return once again to our toy model. In this case, 
important insight may be gained by understanding the specifics of possible thermalization 
processes, as well as more precisely understanding the emergence of quantum mechanics. 

Finally, we point out one more interesting feature of our model. Our scenario required 
the bosonic string, with it's associated cubic instability, in order to dynamically approach the 
double scaling limit. This could indicate that supersymmetry, while important for controlling 
quantum fluctuations in the critical string, has no preferred role in our scenario. Indeed it 
may be the case that apparent instabilities provide critical elements to the dynamics. In this 
context, it may be interesting to consider more general non-critical and non-supersymmetric 
string theories [Zj. 
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